Role of covalent and metallic intercalation on the electronic properties of epitaxial 

graphene on SiC(OOOl) 
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We present an orbital-resolved density functional theory study on the electronic properties of 
hydrogen and lithium intercalated graphene grown on the Si face of SiC. Starting from the (6\/3 x 
GVS)R30° surface reconstruction of the graphene/SiC heterosystem, we find that both H and Li can 
restore the ideal structural characteristics of the two nonequivalent junction parts (i.e. graphene and 
the SiC substrate) when inserted at the interface. However, the chemical/electrostatic interactions 
remain different for the two cases. Hence, H-intercalated epitaxial graphene is subject to a sublattice 
symmetry-breaking electronic interference that perturbs the Dirac point, whereas Li intercalation 
gives rise to a highly n-doped system due to a nonuniform derealization of Li charges. Results 
bring to discussion the role of substrate engineering in epitaxial graphene on SiC. 

PACS numbers: 81.05.ue,73.22.Pr 
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I. INTRODUCTION 

Graphene is a single-atom-thick crystal of sp 2 - 
hybridized carbons that exhibits exceptional electronic 
properties due to its linear dispersion relation around the 
Dirac pointi During the last years the quest for its use in 
fundamentals and applications has boosted the research 
on processes that could allow for a controlled synthesis of 
ordered graphene layers. Within this context, epitaxial 
growth on SiC substrates has emerged as one of the prin- 
cipal technologies for large-scale graphene fabrication 2 -^. 
The process could be better described as a "de-growth" 
one, since SiC crystals are heated up to temperatures 
that allow for an extremely higher sublimation rate of 
the Si atoms with respect to the C ones. The remaining 
C surface atoms recompose to form thin graphite films. 
A proper calibration of the process parameters^ can al- 
low for the formation of single graphene layers directly 
grown on a semi-insulating substrate, i.e. with no need 
to be transfered elsewhere after the end of the thermal 
treatment. 

Epitaxial graphene can grow on both the Si and C 
faces of the (0001) surface of hexagonal SiC polytypes 5 , 
or similarly at the (111) surface of cubic SiC^. Even if 
both theorjs££ and experiment^ argue that C face epi- 
taxial graphene should electrically suffer less interfer- 
ences from the substrate, Si face growth results slower^. 
Hence, a better control can be achieved on the forma- 
tion of single or few-layer structures. In this case a 
first carbon-rich interface layer with a (6-\/3 x 6\/3)-R30 o 
surface reconstruction (often called the buffer layer) 
strongly binds to the substrate and is the precursor 
of the overlying graphene films-J^. Experimental evi- 
dence shows that the presence of this carbon layer has 
a negative influence in the conduction properties of the 
heterosystem with respect to the 5W2-deposited case: 
measurements report an enhanced surface polar phonon 
scattering mechanism^, reduced mean-free paths^- and 
temperature-dependent mobilities^, which are charac- 
teristic of the diffusive transport regime. To this end, 



intercalation techniques with various elements have been 
proposed in the literature^— in order to detach the 
buffer layer from the substrate and minimize the interface 
interaction. A general characteristic of these processes is 
that the intercalated element can form either covalent 
(e.g. H or Ge) or metallic (e.g. Li or Au) bonds with 
the SiC(OOOl) surface and induce a doping level, whose 
origin is often unclear. 

In this article we study the structural and electronic 
properties of H and Li-intercalated epitaxial graphene 
within a comparative approach of functionalization with 
adatoms that interact in a covalent and a metallic way 
with the substrate. As a starting point we recognize 
the problems related to the presence of the buffer layer 
within a density functional theory (DFT) description of 
the (6V3 x 6V3)i?30° reconstruction (Sec. HJ). We then 
proceed with an analysis of the structural and electronic 
symmetry of the intercalated systems (Sec. IIVI and IVl for 
H and Li respectively) , arguing that the resulting chemi- 
cally modified interfaces are different in terms of effective 
doping, bonding and electrostatic interactions. Finally, 
in Sec. |VI]we discuss the versatility of substrate engi- 
neering in tailoring the physical properties of epitaxial 
graphene on SiC. 



II. METHODOLOGY 

We use the DFT SIESTA computational code 18 to per- 
form ab initio calculations, treating the electronic corre- 
lations within the local density approximation (LDA) 19 , 
which has been shown to describe properly the experi- 
mental aspects of the graphene/SiC interfaceiS. More- 
over, taking into account the absence of nonlocal disper- 
sive interactions in the LDA, we also perform a compar- 
ative analysis for the intercalated systems using different 
exchange-correlation functionals: in the spirit of Ref. [2(| 
apart from the LDA we use the Perdew-Burke-Ernzerhof 
implementation 21 of the generalized gradient approxima- 
tion (GGA) and the van der Waals (vdW) functional of 



Dion et a / 22 ' 23 . The epitaxial graphene structures com- 
prise of two bilayers of a (6\/3 x 6\/3)i?30° SiC sub- 
strate (passivated with H at the bottom of the slab), 
over which a single (13 x 13) graphene supercell relaxes 
in order to satisfy lattice commensuration. The aspect 
of the correct surface reconstruction is fundamental in 
the computational modeling of the designated systems. 
Indeed, models based on the simpler (\/3 x \/3)R30 o 
reconstruction 2 ^—, albeit efficient in describing some of 
the basic electronic properties of the SiC /graphene inter- 
face (e.g. n-type doping and presence of the carbon-rich 
layer) , are characterized by: (a) the lack of commensura- 
tion for the graphene and the SiC lattices that leads to 
an 8% stretching of the graphene sheet, and (b) a small 
periodically-repeated (2 x 2) graphene supercell. Reper- 
cussions of these aspects with respect to the results ob- 
tained by the actual (6-\/3 x 6\/3)i?30 o surface symmetry 
are discussed in Sec. |TV] 



We model three different systems, i.e. the 
graphene/SiC interface before and after intercalation 
with H and Li, considering a full and uniform cover- 
age of the substrate by the intercalating element in the 
form of a monolayer. A basis set of double-^ Sankey- 
type valence orbitals has been used for C, Si and H, 
while polarization orbitals have been added in the case 
of Li. Convergence tests on a (\/3 x \/3)i?30 o -based sys- 
tem have shown that the previous set can capture the 
main band structure features obtained by a more ac- 
curate basis that includes polarization orbitals for all 
elements. We exploit the localized character of the 
atomic orbitals in order to calculate the single-orbital 
contributions in the composite electronic band structure: 
for each nonequivalent /c-point sampled along the closed 
r — > M — > K — » r Brillouin-zone path we calculate 
the wave functions 'J'^k corresponding to the n eigen- 
states. Based on the linear combination of atomic or- 
bitals ty n ,]i = X)i (where Ci are the weighting co- 
efficients and 4>i the- pseudoatomic orbitals) we calcu- 
late the square modulus of the total weight of the sin- 
gle orbitals at the n th eigenstate as: w n ^ = J2j l c j| 2 ' 
where the integration runs over those indexes j < i for 
which we want to calculate the contributions in the elec- 
tronic structure. Contrary to the variational approach 
in the calculation of the valence electronic properties, 
the electronic contribution of the ionic cores is stati- 
cally described with norm-conserving Troulier-Martins 
pseudopotentials2£ that have been tested to accurately 
reproduce the band structure of hexagonal SiC polytypes 
and graphene. The minimization of the electron density 
is achieved by sampling the Brillouin zone with a single- 
point Monkhorst-Pack grid for the (6^ x 6^)^30° 
model and a 7 x 7 x 1 grid for the (a/3 x V3)R30° one. A 
mesh cutoff energy of 350 Ryd has been imposed for real- 
space integration, while all structures have been relaxed 
with a force criterion of 0.06 eV/A. 
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FIG. 1. (Color online) (a) Color map topmost view of the 
(6\/3 x 6\/3)-R30 o -reconstructed first carbon-rich layer of an 
epitaxial graphene system showing the vertical positions of 
the C atoms. Here, a gradual yellow (bright in grayscale) to 
red (dark in grayscale) coloring indicates bigger to smaller 
distances from the substrate, (b) Density of states as a func- 
tion of energy for an ideal 4H-SiC crystal (upper) and the 
graphene/SiC(0001) interface (lower). 



III. THE GRAPHENE/SIC(0001) INTERFACE 

As a first step towards understanding the need for in- 
tercalation in epitaxial graphene systems grown on the 
Si face of SiC, we calculate the structural and electronic 
properties of the (6^/3 x 6^3)^30° graphene/SiC inter- 
face. The relaxation of the corresponding supercell gives 
rise to a corrugated carbon layer with a thickness of ~ 1.5 
A. The surface is characterized by height-varying Moire 
patterns that form hexagons with edges of 10-12 A, which 
correspond to the regions of the carbon layer that are 
more distant from the substrate [see the yellow (bright 
in grayscale) areas in Fig. [Ha)]. This structural de- 
formation with respect to the ideal lattice topology is 
also reflected in terms of the chemical bonding, where an 
interplay between sp 2 and sp 3 interface interactions ap- 
pears. Quantitatively, there exists a strong preferential 
disposition towards the sp 3 C-Si covalent bonding, which 
bounds the carbon layer strongly to the substrate. This 
aspect comes in contrast with the respective interface of 
C face epitaxial graphene^. where interface bonding is in- 
herently weaker, and reflects the weakness of the 7r-bond 
in Si with respect to C. The Si atoms of the (0001) sur- 
face that do not covalently bind with the substrate main- 
tain the characteristics of the polar SiC(0001) surface 
with a pronounced inward relaxation with respect to the 
ideal position, while they introduce dangling-bond bands 
within the SiC bandgap 28 . This last feature character- 
izes the electronic properties of the entire graphene/SiC 
heterostructure, where a half-filled peak in the density of 
states at the Fermi level (Ep) appears [Fig. QJb)-lower] . 
Moreover the 2p z orbitals of the carbon epilayer, which 
originate from C atoms that are not bound to the sub- 
strate, give rise to surface states within the SiC bandgap. 
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FIG. 2. (Color online) (a) Geometry of the relaxed H- 
intercalated epitaxial graphene system, (b) Total band struc- 
ture of the H-intercalated system, (c) Contributions of the 
2p z orbitals of the graphene layer in the total band structure, 
(d) Contributions of the substrate orbitals in the total band 
structure. The z-axis is perpendicular to the graphene plain. 



model 


LDA 


GGA 


vdW 


(V3 x ^#30° 


130 meV 


80 meV 


15 meV 


(6^3 x 6^3)i?30° 


8 meV 


5 meV 


5 meV 



TABLE I. Comparison of the energy gap values (E g ) of the El- 
intercalated epitaxial graphene system for the (v3~X \/3)-R30 o 
and the (6^3 x 6^3)i?30° models within the LDA, GGA and 
vdW exchange-correlation functionals. 



The key point that arises from these simulations is that 
the graphene/SiC interface is electrically active. Such 
feature should have an adverse impact on the device 
potential of Si face epitaxial graphene systems, since it 
could be directly related to the manifestation of leakage 
interface currents during conduction. 

IV. HYDROGEN INTERCALATION 

Si-face grown epitaxial graphene can be used as a ref- 
erence system for the study of a number of fundamental 
physical properties 2 ^ without the need of further process- 
ing. However, the electrical use of such material within 
a device concept dictates the need of an interface pas- 
sivation with functional adatoms through a post-growth 
chemical process. The simplest example of intercalation 
from a monovalent element is that of hydrogen, as pro- 
posed by Riedl et al. 13 and followed by other s 12 ' 30 . Here, 
starting from the (6^/3 x 6-\/3)-R30 o reconstruction of the 
carbon-rich layer, we model H intercalation by assign- 
ing a single H atom on top of each Si atom of the sub- 
strate's surface. This geometry is experimentally con- 
firmed by infrared absorption spectroscopy-^ 2 -. Upon re- 
laxation, H atoms covalently bind with the interface Si 
atoms and restore the ideal sp 3 hybridization of the sub- 
strate [Fig. (Ha)]. As a result, the structurally disordered 



carbon layer turns into an almost detached flat graphene 
sheet with a distance of 2.39 A from the H interface layer 
and 3.92 A from the substrate Si atoms. The electronic 
band structure of this system shows the (13 x 13) folded 
graphene bands and confirms the minimal interaction 
with the substrate (Fig. HJb)). However a small pertur- 
bation appears at the Dirac point that opens a bandgap 
of few millielectronvolts, in analogy with a similar feature 
calculated for the interaction between graphene and the 
(0001) surface of SiO^. We find this perturbation ro- 
bust for small horizontal shifts of the graphene layer with 
respect to the substrate. Even if this gap is extremely 
small, well beyond fc^T at room temperature and difficult 
to reveal in experiments, its physical origin is of a partic- 
ular interest. A careful analysis of the electronic Hamil- 
tonian of the composite system excludes an inter-valley 
contribution on the electronic perturbation. Instead, the 
origin of this interference can be traced at the breaking of 
the sublattice symmetry 3 - 2 , due to the unequal coupling 
between the two interpenetrating triangular lattices of 
the graphene sheet and the localized Si-H dipoles of the 
substrate, which have non-equivalent positions with re- 
spect to the graphene atoms in the supercell. For the 
smaller (V3 x y/3)R30° model, this symmetry-breaking 
effect acquires an extreme form due to the fact that the 
graphene sheet is defined by only eight periodically re- 
produced C atoms, and manifests in the form of wider 
band gaps in the electronic band structure (see Table IJ). 
In this case, the choice of the exchange-correlation func- 
tional becomes critical and both the LDA and the GGA 
overestimate the graphene-substrate interaction. Con- 
trary, results obtained by the vdW functional indicate 
that the inclusion of nonlocal interactions due to fluctu- 
ating dipoles reduces the charge correlations at the inter- 
face of H-intercalated epitaxial graphene systems. When 
considering the (6-\/3 x 6-\/3)i?30 o model with the 338 
graphene atoms in the supercell, the sublattice differ- 
ences get almost averaged out and all functionals con- 
verge in a minimal Dirac-point perturbation, which is 
consistent with the experiment. It should be noted here 
that this effect is inherent with nonuniform electrostatic 
potentials induced by substrates with atomic arrange- 
ments that are unequal to that of graphen o 10 ' 31 . How- 
ever, its quantitative manifestation in the H-passivated 
SiC(0001) obtained by the (6\/3 x 6^)^30° cell, 

is negligible. A second issue worth mentioning for the 
(y/3 x %/3)i?30° case is the shift of the graphene eigen- 
states with respect to the substrate levels due to ficti- 
tious strain, which results in a positioning of the Dirac 
cone closer to the valence band of SiC with respect to the 
bigger supercell. Nonetheless, all models converge in the 
designation of charge neutrality at the Dirac point. This 
aspect can be twofold commented: on one hand the ab- 
sence of an effective doping in this system shows that the 
H-passivated SiC (0001) surface looses the characteristics 
of the free polar SiC (0001) surface and does not introduce 
Fermi-level-pinning midgap states in the electronic struc- 
ture. On the other hand, the simulation of an extremely 
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FIG. 3. (Color online) (a) Side view of the relaxed Li- 
intercalated epitaxial graphene system, (b) Top view of the 
Li-intercalated system, (c) Top view showing only the Li 
atoms that occupy the hollow positions of the first SiC bi- 
layer. 



thin SiC slab neglects the formation of intrinsic dipoles 
formed at the (0001) direction of hexagonal SiC poly- 
types and may undermine their doping contributions. In 
this sense, the results presented in this paragraph should 
be more suitable for cubic, rather than hexagonal SiC. 



V. LITHIUM INTERCALATION 

The main problem with the use of atomic hydrogen 
in the chemical detachment of the graphene/SiC inter- 
face is that it can also functionalize the graphene layer 
itself^, giving rise to sp 3 -type defects. Such defects are 
at the origin of transport pseudogaps arround the Dirac 
point during electrical conduction and can substantially 
increase the electronic localization and, hence, the elec- 
trical disorder—. A possible solution towards this di- 
rection could be the intercalation with elements that do 
not disturb the planarity and the sp 2 character of the 
graphene sheet. Lithium fulfills this prerequisite 3 -^. Ex- 
perimentally Li-intercalation has been proposed by Vi- 
rojanadara et al*^ who showed that Li, as a highly re- 
active lightweight monovalent metal, can penetrate the 
strongly-bound carbon-rich layer and position on top of 
the SiC substrate. Motivated by such experiment, we 
model Li-intercalated epitaxial graphene systems consid- 
ering an equivalent to H complete coverage of the SiC 
surface. Structurally, also in this case the carbon inter- 
face layer relaxes into a flat graphene position with a 
vertical distance of 2.37 A from the Li adlayer and 4.41 
A from the Si atoms of the substrate [Fig. [3Ja)] - How- 
ever, as a first difference with respect to the H case, we 
find that the minimum-energy configuration for the re- 
laxation of Li atoms is at the hollow positions of the first 
substrate bilayer [Fig. 02[b-c)]. The reason of this sym- 
metry originates from the bonding interactions between 
the Li adlayer and the substrate. We find an important 
charge transfer from the Li atoms to a region between 
Li and the Si atoms of the substrate's surface, localized 



on the top of the Si atomic positions. Such configura- 
tion guarantees the pure sp 3 hybridization of the sub- 
strate, while it leaves a cation at the location of the Li 
atom. The unequal chemisorption of H and Li on the 
graphene/SiC interface could be also reflected in a differ- 
ent sub-graphene inter-diffusion mechanism during the 
intercalation process. 

The electronic band structure obtained by this com- 
posite system shows highly-concentrated bands through- 
out the hole region of the spectrum that are difficult to 
interpret. In this sense an orbital-resolved analysis is nec- 
essary for the determination of the local contributions 
in the total electronic structure. Considering therefore 
the contribution of the 2p z orbitals that originate from 
the graphene layer, we obtain a band structure that is 
equal to that of ideal graphene for a large energy range 
around the Dirac point [Fig. Ufa)]. This result denotes 
that the detachment from the substrate also in the case 
of Li intercalation is almost complete, while unlike the 
H-intercalated case we do not obtain any interference at 
the Dirac point due to a better screening of the substrate 
from the Li ionic cores. The important aspect though is 
that the graphene sheet is highly electron-doped (as also 
calculated for a low-coverage Au intercalation^), with 
the Dirac cone lying ~ 1.3 eV below the Fermi level of 
the heterosystem. This high effective doping is due to the 
electrostatic interactions between the graphene layer and 
the Li cations, which gives rise to a shift of the graphene 
bands within the valence band region of the SiC sub- 
strate [Fig. IHc)], making the whole system electrically 
active at the vicinity of the Dirac point. The analysis of 
the electronic Hamiltonian of the composite system here 
shows an important decrease of the on-site energies of the 
C atoms with respect to the case of H-intercalation. This 
aspect is important if we consider that the actual interca- 
lation process may not be uniform throughout the entire 
breadth of the graphene layer. In this case we can ex- 
pect an Anderson-type disorder arising from the energetic 
mismatch between intercalated and non-intercalated re- 
gions. A further proof of the differences between the H 
and Li-intercalated systems results from the calculation 
of the binding energies of the graphene layer on the pas- 
sivated substrate: here the LDA values are 0.06 eV and 
0.108 eV per C atom respectively. Finally, it is interest- 
ing to note that the quasi-flat bands originating from the 
Li orbitals [Fig. Hfb)] maintain a higher-energy position 
with respect to the Dirac cone and, hence, should not 
interfere during electrical conduction around Dirac-point 
energies. 



VI. DISCUSSION 

The extrapolation of device-requested characteristics, 
like high electron mobilities and densities, from epi- 
taxial graphene systems grown on the Si face of SiC 
makes necessary the passivation of the interface between 
graphene and SiC. Towards this objective, intercalation 



FIG. 4. (Color online) Orbital-resolved band structure for the Li-intercalated epitaxial graphene system, showing the contri- 
butions of the (a) 2p z orbitals originating from the graphene C atoms, (b) all orbitals originating from Li atoms and (c) all 
orbitals originating from the substrate C/Si atoms. The z-axis is perpendicular to the graphene plain. 



techniques with various functional elements have been 
proposed in the literature. In this article we have stud- 
ied within a DFT framework the electronic structure 
properties of the non-intercalated as well as interca- 
lated interfaces between graphene and the Si face of SiC, 
starting from the actual (6-\/3 x 6\/3)-R30 o reconstruc- 
tion of the graphene/SiC heterosystem. We identified 
some issues related to the graphene/SiC coupling, like 
the presence of midgap interface states within the SiC 
bandgap, which could give rise to plausible faults in the 
expected electrical behavior. We thereon studied two 
different intercalated systems with H and Li adatoms 
respectively, and saw that although in both cases the 
strongly-bound carbon-rich layer gets detached from the 
substrate, different electronic properties distinguish the 
two structures: in the case of H intercalation we observed 
a minimal perturbation at the Dirac point as a conse- 
quence of substrate-induced interference from localized 
Si-H dipoles, while in the case of Li intercalation we ob- 
tained a highly n-doped ideal graphene band structure 
with the position of the Dirac cone within the valence 
band of the SiC substrate. A critical comparison between 



these two examples unfolds the potentiality of substrate 
engineering in epitaxial graphene on SiC. 

It can be argued that research on intercalation pro- 
cesses for the use of epitaxial graphene in applications 
should focus on: (a) a full and uniform substrate cov- 
erage, (b) an absence of functionalizcd areas on the 
graphene layer, and (c) an absence of non-graphene states 
around the Dirac point. Moreover, the level of doping can 
be appropriately used in order to engineer the active re- 
gions for device components. Within this context, both 
covalent and metallic intercalation mechanisms could be 
useful for future epitaxial graphene device integration. 
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